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Signature kinematic variables and characteristic concentrations in phase space of double parton 
scattering are discussed. These properties should allow the double-parton contribution to 
pp — > b b jet jet X at Large Hadron Collider energies to be distinguished from the usual 
single parton scattering contribution. A methodology is suggested to measure the size of the 
double-parton cross section. 

1 Introduction 

Double parton scattering (DPS) means that two short-distance hard-scattering subprocesses 
occur in a given hadronic interaction, with two initial partons being active from each of the 
incident protons in a collision at the Large Hadron Collider (LHC). The concept is shown for 
illustrative purposes in the left diagram of Fig. [TJ It may be contrasted with conventional single 
parton scattering (SPS) in the right diagram, in which one short-distance subprocess occurs, 
with one parton active from each initial hadron. Both contribute to the same 4 parton final 
state. Processes such as sketched in the left diagram of Fig. [T] are included in descriptions of the 
underlying event in some Monte Carlo codes. Our interest is to investigate whether this second 
hard process can be shown to be present in LHC data, as a perturbatively calculable hard part 
of the underlying event. 

Studies of double parton scattering have a long history theoretically, with many references to 
prior work listed in our paper^, and there is evidence in data^. A greater role for double-parton 
processes may be expected at the LHC where higher luminosities are anticipated along with the 
higher collision energies. A large contribution from double parton scattering could result in a 
larger than otherwise predicted rate for multi-jet production, and produce relevant backgrounds 
in searches for signals of new phenomena. The high energy of the LHC also provides an increased 
dynamic range of available phase space for detailed investigations of DPS. 

Our aims^ are to address whether double parton scattering can be shown to exist as a 
discernible contribution in well defined and accessible final states, and to establish the charac- 
teristic features that allow its measurement. We show that double parton scattering produces 
an enhancement of events in regions of phase space in which the contribution from single parton 
scattering is relatively small. If such enhancements are observed experimentally, with the kine- 
matic dependence we predict, then we will have a direct empirical means to measure the size 
of the double parton contribution. In addition to its role in general LHC phenomenology, this 
measurement will have an impact on the development of partonic models of hadrons, since the 
effective cross section for double parton scattering measures the size in impact parameter space 
of the incident hadron's partonic hard core. 



Figure 1: (diagram on the left) Sketch of a double-parton process in which the active partons are i and k from one 
proton and j and I from the second proton. The two hard scattering subprocess are A(i j —> ab) and B(k I — > c d). 
(diagram on the right) Sketch of a single-parton process in which the active partons are i from one proton and j 
from the second proton. The hard scattering subprocess is A(i j — > a b c d). 



From the perspective of sensible rates and experimental tagging, a good process to examine 
should be the 4 parton final state in which there are 2 hadronic jets plus a b quark and a b 
antiquark, viz. b b j\ j%. If the final state arises from double parton scattering, then it is 
plausible that one subprocess produces the b b system and another subprocess produces the two 
jets. There are, of course, many single parton scattering (2 to 4 parton) subprocesses that can 
result in the b b j\ j 2 final state, and we identify kinematic distributions that show notable 
separations of the two contributions. 

The state-of-the-art of calculations of single parton scattering is well developed whereas the 
phenomenology of double parton scattering is less advanced. For pp — > bbjij 2 X, assuming that 
the two subprocesses A{i j — > b b) and B(k I — > j\ j 2 ) in Fig. [I] are weakly correlated, and that 
kinematic and dynamic correlations between the two partons from each hadron may be safely 
neglected, we employ the common heuristic expression for the DPS differential cross section 

a OPS, _^iX- ■ Y \ da SPS (pp ->• bbX)da SPS {pp -»• jij 2 X) 

da °{pp bb]ij 2 X) = . (1) 
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The numerator is a product of single parton scattering cross sections. In the denominator, there 
is a term <T e fj with the dimensions of a cross section. Given that one hard scatter has taken place, 
c c ff measures the effective probability for a second hard scatter. Collider data^l yield values in 
the range a e g ~ 12 mb. We use this value for the estimates we make, but we emphasize that 
the goal should be determine its value in experiments at LHC energies. 

The details of our calculation of the double parton and the single parton contributions to 
p p — > bb j\ J2 X are found in our paper^. We perform full event simulations at the parton level 
and apply a series of cuts to emulate experimental analyses. We also treat the double parton 
and the single parton contributions to 4 jet production, again finding that good separation is 
possible despite the combinatorial uncertainty in the pairing of jets^. 

2 Distinguishing variables 

Correlations in the final state are predicted to be quite different between the double parton 
and the single parton contributions. For example, we examine the distribution of events as 
function of the angle $ between the planes defined by the bb system and by the jj system. 
If the two scattering processes ij — > bb and kl — > jj which produce the DPS final state are 
truly independent, one would expect to see a flat distribution in the angle By contrast, 




Figure 2: (left panel) Event rate as a function of the angle between the two planes denned by the bb and jj 
systems. In SPS events, there is a correlation among the planes which is absent for DPS events, (right panel) 
The transverse momentum pr distribution of the leading jet in jjbb, either a b jet or a light jet j. 



many diagrams, including some with non-trivial spin correlations, contribute to the 2 parton to 
4 parton final state in SPS ij — > 66jetjet, and one would expect some correlation between the 
two planes. In the left panel of Fig. [21 we display the number of events as a function of the 
angle between the two planes. There is an evident correlation between the two planes in SPS, 
while the distribution is flat in DPS, consistent with the expectation that the two planes are 
uncorrelated. 

Another dynamic difference between DPS and SPS is the behavior of event rates as a function 
of transverse momentum. As an example of this, in the right panel of Fig. [21 we show the 
transverse momentum distribution for the leading jet (either a b or light j) for both DPS and SPS. 
SPS produces a relatively hard spectrum, associated with the presence of several propagators in 
the hard-scattering matrix element. On the other hand, DPS produces a much softer spectrum 
which (up to issues of normalization in the form of <7 e fj) can dominate at small values of transverse 
momentum. For the value of a e g and the cuts that we use, SPS tends to dominate over the full 
range of transverse momentum considered. The cross-over between the two contributions to the 
total event rate is ~ 30 GeV for the acceptance cuts considered. A smaller (larger) value of <7 e ff 
would move the cross-over to a larger (smaller) value of the transverse momentum of the leading 
jet. 

Although interesting, the two distributions in Fig. [2] would not allow the two components, 
DPS and SPS, to be separated. We turn next to the search for variables that could allow a 
clear separation of the contributions. At lowest order for a 2 — > 2 process, the vector sum 
of the transverse momenta of the final state pair vanishes, although in reality, radiation and 
momentum mismeasurement smear the expected peak near zero. Nevertheless, the DPS events 
are expected to show a reasonably well-balanced distribution in the transverse momenta of the 
jet pairs. To encapsulate this expectation for both light jet pairs and 6-tagged pairs, we use the 
variable^: 

q , J_ If \Mh,b2)\ \ 2 ( Mii,j 2 )| \ 2 ( ) 

PT V2 V \\PT(h)\ + MW \\Pr(h)\ + \pt{32)\) " 1 ' 

Here pr(bi,b2) is the vector sum of the transverse momenta of the two final state b jets, and 
Pr(ji, 32) is the vector sum of the transverse momenta of the two (non b) jets. We expect pt ~ 
for both of these vector sums. 

The distribution in S' p is shown in Fig. EJa). As expected, the DPS events are peaked near 
S' pT ~ and are well-separated from the total sample. The SPS events, on the other hand, tend 
to be broadly distributed and show a peak near S' ~ 1. The peak near 1 is related to the fact 




Figure 3: (a) Distribution of events in S' PT for the DPS and SPS samples. Due to the back-to-back nature of 
the 2 — > 2 events in DPS scattering, the transverse momenta of the jet pair and of the fr-tagged jet pair are 
small, resulting in a small value of S' PT . (b) The variable S 1 ^ for DPS and SPS+DPS events provides a stronger 
separation of the underlying DPS events from the total sample when compared to A(j> for any pair. 

that a significant number of the SPS bb or jj pairs arise from gluon splitting which yields a large 
Pt imbalance and, thus, larger values of S' . 

The azimuthal angle between pairs of jets is another variable that can represent the roughly 
back-to-back hard-scattering topology of the DPS events. We expect the azimuthal angle be- 
tween the pairs of jets corresponding to each hard scattering event to be strongly peaked near 
A^jj ~ ~ 7T. Real radiation of an additional jet, where the extra jet is missed because 

it fails the threshold or acceptance cuts, allows smaller values of Ac^. There is a clear peak 
near Acftjj = it for DPS events, while the events are more broadly distributed in SPS events^. 
A secondary peak near small A.(pjj arises from gluon splitting which typically produces nearly 
collinear jets. As in the case of the S' pT variable, the separation of DPS events from SPS events 
becomes more pronounced if information is used from both the bb and jj systems. One dis- 
tribution built from a combination of the azimuthal angle separations of both jj and bb pairs 

isEl 

S4, = 4 V / A0(6i,6 2 ) 2 + A0(j 1 ,j 2 )2. (3) 

The SPS events are broadly distributed across the allowed range of S^, shown in Fig. 0(b). 
However, the combined information from both the bb and jj systems shows that the DPS events 
produce a sharp and substantial peak near ~ ir which is well-separated from the total 
sample. The narrow peaks near Acftjj = ir and near Sj, = 1 are smeared somewhat once soft 
QCD radiation and other higher-order terms are included in the calculation. 

In our simulations, the variable S' pT appears to be a more effective discriminator than S^. 
However, given the leading order nature of our calculation and the absence of smearing associated 
with initial state soft radiation, this picture is subject to change, and a variable such as (or 
some other variable) may offer a clearer signal of DPS at the LHC. Realistically, it would be 
valuable to study both distributions, once LHC data are available, in order to determine which 
is more instructive. 

The evidence in one-dimensional distributions for distinct regions of DPS dominance prompts 
the search for greater discrimination in a two dimensional distribution of one variable against 
another. One scatter plot with interesting features is displayed in Fig. [H The DPS events are 
seen to be clustered near SL = and are uniformly distributed in <£. The SPS events peak 
toward S pT = 1 and show a roughly sin<i> character. While already evident in one-dimensional 
projections, these two features are more apparent in the scatter plot Fig.[H Moreover, the scatter 
plot shows a valley of relatively low density between 5' ~ 0.1 and ~ 0.4. In an experimental 
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Figure 4: Two-dimensional distribution of events in the inter-plane angle $ and the scaled transverse momentum 

variable S' PT for the DPS and SPS samples. 

one-dimensional distribution, one would see the sum of the DPS and SPS contributions. If 
structure is seen in data similar to that shown in the scatter plot Fig. HJ one could make a cut 
at S' pT < 0.1 or 0.2 and verify whether the experimental distribution in <I> is flat as expected for 
DPS events. 

3 Strategy and Further Work 

The clear separation of DPS from SPS events in Fig. [J] suggests a methodology for the study of 
DPS. One can begin with a clean process such as pp — > bbjij2X and examine the distribution 
of events in the plane defined by S' pT and $. We expect to see a concentration of events near 
S' pT = that is uniformly distributed in <£. These are the DPS events. Assuming that a valley 
of low density is observed between S' pT ~ 0.1 and ~ 0.4, one can make a cut there that produces 
an enhanced DPS sample. Relative to the overall sample, this enhanced sample should show a 
more rapid decrease of the cross section as a function of the transverse momentum of the leading 
jet, and the DPS enhanced sample can be used to measure a e g. A similar examination of other 
final states, such as 4 jet production, will answer whether the extracted values of <7 e fr are roughly 
the same. Theoretical and experimental studies of other processes can follow, such as bbti, Wjj, 
and Hjj. 

On the phenomenological front, next-to-leading order (NLO) expressions should be included 
for both the SPS and DPS contributions. The NLO effects are expected to change normalizations 
and, more importantly, the distributions in phase space. The sharp peaks near — 7r in 
Fig. [3](b) and S' = in Fig. [3]^a) will be broader and likely displaced somewhat. The weak 

correlation between subprocesses assumed in Eq. (pQ) cannot be strictly true^. With a large 
enough data sample at the LHC one could investigate the extent to which correlations play a 
significant role. 

Finally, it would be good to examine the theoretical underpinnings of Eq. ([1]) and, in the 



process, gain better insight into the significance of a c g. A firm basis is desirable for Eq. (pQ) 
starting from the formal expression for the differential cross section in terms of the absolute 
square of the full matrix element integrated over phase space: 

da( PP -> bbj lj2 X) = j- 1 M(pp -> bbjij2X)\ 2 dPS bihnX . (4) 

The amplitude M(pp — > bbjij 2 X) should include a sum of amplitudes for 2-parton collisions 
(one active from each incident hadron, i.e., 2—7-4); 3-parton collisions (two active from one 
hadron and one active from the other); and 4-parton collisions (two active from each hadron 
or three from one and one from the other), and so forth that all yield the same 4 parton final 
state. There will be contributions to the final state from the squares of individual amplitudes 
as well as interference terms. Specializing to 4 — > 4, the DPS case, one would start from a 
4-parton — > 4-parton hard part. Not evident at this time is how the four-parton matrix element 
can be reduced to a product of two matrix elements for the single parton scatterings, needed for 
Eq. (pQ). The demonstration of clear DPS signals in LHC data would be an important stimulus 
for further theoretical studies. 
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